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line). Such scaling matches well with the theo- 
retical prediction of optical Stark shift 5(^coo) = 
2S'-Ep /hQ., where S is the optical Stark shift con- 
stant related to the transition dipole moment 
and Ep is the electric field of the pump pulse 
{25, 28). From the experimental data, we can de- 
termine an optical Stark shift constant S ~ 45 
for A-exciton in monolayer WSe2, which is of 
similar magnitude to that for exciton transition 
in semiconductor quantum wells (29). 

The valley-selective optical Stark shift breaks 
the degeneracy of valley exciton transitions in 
monolayer WSe2 and defines an effective valley 
pseudo-magnetic field. In our experiment, the 
photo-induced energy splitting between K and 
K' exciton transitions can be as large as 10 meV 
(Fig. 4 and supplementary text). The correspond- 
ing pseudo-magnetic field ^eff for valley exci- 
tons can be 

estimated with ^eff = , where |Ib is the 

Bohr magneton and gex is the effective ^--factor 
for valley exciton transitions in WSe2. The effective 
exciton ^--factor g^^ combines contributions from 
both electrons and holes and has a theoretically 
predicted value of -1.5 for WSe2 (supplementary 
text). Using this g'-factor, we estimate a pseudo- 
magnetic field ^eff as high as 60 T for a 10-meV 
splitting of valley exciton transitions. A real mag- 
netic field of this magnitude is difficult to achieve 
even with superconducting magnets, but such a 
pseudo-magnetic field for MX2 valley excitons 
can be produced conveniently and with femto- 
second temporal control by using light pulses. 

It has been reported recently that excitons in 
different valleys in monolayer WSe2, resonantly 
excited by linear polarization light, can maintain 
their phase coherence over extended time (7). 
The valley-dependent optical Stark effect offers 
a convenient and ultrafast way to lift the valley 
degeneracy and introduce a controlled phase dif- 
ference between the two valley states, therefore 
enabling coherent rotation of resonantly excited 
valley polarizations with high fidelity (fig. S2). In 
analogy with spintronics, such coherent manip- 
ulation of valley polarization can open up fasci- 
nating opportunities for valleytronics. 
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HEAVY FERMIONS 



The recent development of topological insu- 
lators is a triumph of single electron band 
theory (1-8). Kondo insulators can be used 
to explore whether similar exotic states of 
matter can arise in the presence of strong 
electronic interactions. In these strongly correla- 
ted heavy-fermion systems (9), the hybridization 
between itinerant electrons and localized orbi- 
tals opens a gap and makes the material insula- 
ting. Once the sample temperature is cold enough, 
the electronic structure can be mapped to a state 
that resembles a normal topological insulator (TI) 
(10), with a bulk insulating state and a conductive 
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surface state. In samarium hexaboride (SmBg), 
the existence of the surface state has been sug- 
gested by recent experimental observations of the 
surface conductance as well as a mapping of the 
hybridization gap (11-13). We report the obser- 
vation of quantum oscillations in Kondo insu- 
lator SmBg using torque magnetometry via the 
de Haas-van Alphen (dHvA) effect. The observed 
Fermi surfaces are shown to be two-dimensional 
(2D) and arise from the crystalline (101) and 
(100) surfaces. 

One major difference between SmBg and the 
conventional topological insulators is the crystal 
structure, which for SmBg is simple cubic (Fig. 
lA). SmBg single crystals were grown by conven- 
tional flux methods. Each sample was etched with 
acid to remove the leftover flux. Figure IB shows 
a photo of a piece of SmBg single crystal. Beside a 
flat (001) surface, there are four (101) planes. 

We use torque magnetometry to resolve the 
Landau Level quantization and the resulting 
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In the Kondo insulator samarium hexaboride (SmBe), strong correlation and band 
hybridization lead to an insulating gap and a diverging resistance at low temperature. 
The resistance divergence ends at about 3 kelvin, a behavior that may arise from surface 
conductance. We used torque magnetometry to resolve the Fermi surface topology in 
this material. The observed oscillation patterns reveal two Fermi surfaces on the (100) 
surface plane and one Fermi surface on the (101) surface plane. The measured Fermi 
surface cross sections scale as the inverse cosine function of the magnetic field tilt 
angles, which demonstrates the two-dimensional nature of the conducting electronic 
states of SmBe. 
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quantum oscillations in magnetization. Electro- 
nic transport measurements have been used to 
detect quantum oscillations in conventional to- 
pological materials (14-22). In contrast, magne- 
tization is simply the derivative of the magnetic 
free energy G with respect to magnetic field 
H. Therefore, torque magnetometry probes the 
oscillations in the free energy and is sensitive to 
Fermi surface (FS) topology of both 3D and 2D 
electronic systems (23-27). Torque magnetom- 
etry measures directly the anisotropy of the 
magnetic susceptibility of a sample (25). With the 
tilted magnetic field H confined to the d - c 
plane, the torque i of a paramagnet is shown as 
follows: 

= |iQAxi^'^sin(t)COS(t) (1) 

where |io is the vacuum permeability, (\> is the tilt 
angle ofH away from the crystalline a - c axis, 
and Ax = Xa - Xc is the magnetic susceptibility 
anisotropy. Therefore, any change of the FS to- 
pology caused by the Landau Level quantization 
is revealed by torque magnetometry. 



In our experimental setup, an SmBg single crys- 
tal is glued to the tip of a thin brass cantilever 
(Fig. IC, inset). The magnetic torque i is mea- 
sured by tracking the capacitance change be- 
tween the cantilever and a gold film underneath. 
An example torque i versus magnetic field \ioH 
curve (Fig. IC), taken at temperature T = 0.3 K 
and (j) ~ 44°, is quadratic overall, refiecting the 
linear H dependence of the sample magneti- 
zation and the paramagnetic nature of SmBg. 
Large oscillations and small wiggles start to 
appear as the magnetic field goes beyond 5 T; 
for [LqH > 10 T, the fast oscillation patterns 
dominate. 

The first question that arises is whether it is 
possible for the thin atomic layer of the surface 
state to be responsible for the observed mag- 
nitude of the magnetic torque and magnetiza- 
tion. The magnitude of the observed magnetic 
torque was recorded by the relative capacitance 
change, then converted to absolute values using 
the calibrated spring constant of the cantilevers 
(see supplementary materials). Figure ID shows 
the field dependence of the effective magnetic 
moment M= in a magnetic field H as high 
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Fig. 1. Experimental setup. (A) The crystal structure of Kondo insulator SmBe. The Be cluster and rare 
earth element Sm form a simple cubic structure. (B) Photograph of an SmBe single crystal. Guide lines 
are drawn to label the crystalline plane of each surface. (C) The field dependence of the magnetic torque 
X in SmBe. The data were taken at temperature T = 300 mK. Inset: sketch of the measurement setup. 
The sample stage is rotated to tilt magnetic field H in the crystalline a - c plane. The magnetic field was 
applied to the sample with a tilt angle <\) -44° relative to the crystalline c axis. The spring constant of the 
cantilever was determined using the sample weight. (D) The field dependence of the magnetic moment 
M, taken at 0.3 K and under a magnetic field as high as 45 T. The low-field linear M-H dependence 
reflects the paramagnetism in SmBe. The vertical gray scale bar shows the M value corresponding to 
0.5 |iB per surface unit cell. 



as 45 T. The maximum oscillatory M is around 
3 X 10"^^ A-m^. Using the total area of all the 
surfaces of the SmBg sample, including the 
dominating top and bottom (100) surfaces and 
other small edge surfaces, we can estimate that 
there are 7.5 x 10^^ unit cells on the sample 
surfaces. By normalizing the magnetic moment 
M to the number of surface unit cells, we find 
the maximum change of oscillatory Mis around 
AM ~ 0.4 |Ib per surface unit cell, where |Ib is the 
Bohr magneton. The gray scale bar in Fig. ID 
marks 0.5 |Ib per surface unit cell for comparison. 

Our measured oscillatory M is consistent 
with the surface carrier density measured by 
the Fermi pocket size. As we will demonstrate 
later, the measured pocket size (300 T ~ 400 T) 
leads to a total carrier density of ~1 x lO^Vcm^, 
which is -0.17 electrons per unit cell. The the- 
oretical value of the oscillatory M is 2 |Ib per 
electron in 2D electronic systems, confirmed 
experimentally in (26). The carrier density im- 
plies a magnetization oscillation of 0.34 |Ib per 
surface unit cell. Our observed maximum AM ~ 
0.4 |Ib per surface unit cell is consistent with this 
simple estimate. 

Figure 2 presents the oscillation pattern of the 
magnetic torque. The oscillatory torque losc is 
defined by subtracting a quadratic background 
from the torque i. losc is periodic in l/|io^ (Fig. 
2A), reflecting the quantization of the Landau 
Levels. The magnetic torque is measured in 77 up 
to 45 T. For metals, the oscillation frequency F is 
determined by the cross section area A of the 
Fermi surface (23)'. 



where h is the reduced Planck constant, and e 
is the electrical charge. The fast Fourier trans- 
formation (FFT) of the oscillatory toscj Plea- 
sured at 0.3 K and cj) - 32° (Fig. 2B) shows a 
number of peaks, including a small pocket a at 
i^" - 35 T, a larger pocket ^^\F^ ~ 300 T, and 
the largest pocket y at i^^ ~ 450 T. In the higher- 
frequency range, we observed a series of new 
features: the second harmonic of the p pocket 
at 2p ~ 600 T and the third harmonic of the p 
pocket at 3P ~ 900 T. Most notably, the other 
peak (marked as appears at -520 T and 
arises from the oscillation due to the p pocket 
in a neighboring (101) surface plane. This deter- 
mination is based on the angular dependence of 
the oscillation frequencies, as discussed in detail 
below and in Fig. 3. 

The electronic properties of these three Fermi 
surfaces are revealed by tracking the temper- 
ature and field dependence of the oscillatory 
torque losc (Fig- 2C). Even though an insulat- 
ing gap is observed in the temperature depen- 
dence of the conductivity of SmBg (11-13), the 
temperature dependence of Tqsc is very much 
like that of a normal metal. In metals, the os- 
cillating magnetic torque is well described by 
the Lifshitz-Kosevich (LK) formula (23). The 
temperature and magnetic field dependences 
of the oscillation amplitude are determined 
by the product of the thermal damping factor 
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i?T and the Dingle damping factor R^,, de- 
fined as 



Rt = aTm*/Bsmh{aTm*/B) 
Rd = exp{-aTDm*/B) 



(3) 



where the effective mass m = m*mg and the 
Dingle temperature = h/2nkBTs ■ is the 
scattering rate, is the Boltzmann Constant, 
TRe is the bare electron mass, B = [LqH is the 
magnetic flux density, and a = 2n^k^meleh 
-14.69 T/K {23). 

Fitting the temperature dependence of the 
normalized oscillation amplitudes to the ther- 
mal damping factor R^ (Fig. 2C) yields m = 
0.119me for Fermi surface a, m = 0.129me for 
Fermi surface p, and m = 0.192me for Fermi sur- 
face y. Figure 2D displays the field dependence 
of the oscillation amplitude, normalized by the 
thermal damping factor i?x- Fitting the curves 
to i?D yields Tj) -15.9 K for Fermi surface a, Tjy 
-18.6 K for Fermi surface % and Td - 29.5 Kfor 
Fermi surface y. Based on the extracted oscilla- 
tion frequencies, effective masses, and Dingle 
temperatures, we are able to characterize the 
observed Fermi surfaces in SmBg (Table 1). 

The long mean free path I measured for all three 
pockets (Table 1) supports the idea of intrinsic 
metallic surface states. The mean free path is 
two orders of magnitude larger than the crystal 
lattice constant of -0.4 nm. Such a long mean 
free path is incompatible with the surface/bulk 
conductance arising from hopping between 
impurities. 

The low effective mass that we measured 
(Table 1) is quite surprising, because most of the 
theoretical work on the topological Kondo in- 
sulator predicts a heavy mass of the surface 



states (28). The observed mass in SmBg is also 
much smaller than that of the divalent hexa- 
boride compounds such as La-doped CaBg and 
EuBg (29, 30). Pockets with heavier masses may 
still exist but could not be observed in our ex- 
perimental conditions because of the thermal 
and Dingle dampings. Measurements in higher 
magnetic field and in colder temperatures would 
help resolve such extra Fermi pockets. 

The most important feature of the observed 
quantum oscillations is the two-dimensionality 
(2D). For a 2D planar metallic system, the mag- 
netic field projected to the normal axis of the 
plane determines the Landau Level quantizations. 
Thus, the oscillation frequency versus tilt angle 
curve generally follows the inverse of a sinusoi- 
dal function (15-17). We rotate the whole canti- 
lever setup to track how the oscillation frequency 
F changes as a function of the tilt angle 4) in a 
broad range of more than 180°. Using the com- 
plete FPT plots (figs. S3 to S5), we obtain the an- 
gular dependences of the oscillation frequencies 
of all the Fermi pockets on the (100) and (101) 
family surfaces, as displayed in Fig. 3, A to C. 

For the dominating oscillatory pattern from 
pocket p, the oscillation frequency not only 
displays a large angular dispersion but also 
closely tracks the 2D angular dependence from 
the (101) surface families. There are four branches 
of the - 4) patterns, as a result of the fourfold 
crystalline symmetry of the SmBg cubic struc- 
ture. At each ^, there are two appearances of 
the p family. For example, at cj) - 30° there is one 
F^ from the (101) plane and another one from 
the neighboring (101) plane (the latter contri- 
bution is marked as p' in Fig. 2B). Most notably, all 
the points track the solid lines of the function 
i^oVcosrt + 45°),FoVcos((^-45°),FoVcos((t)-135°), 



and F^/cos((t)-225°). We observe i^"^ at frequen- 
cies as high as 900 T— a more than 200% increase 
from the minimum value of 286 T. Such large 
divergence and the close tracking of the inverse 
cosine dependence strongly support the 2D na- 
ture of the observed p pocket on the (101) surface 
plane families. 

In contrast, the angular dependence of the 
oscillation frequency of the a pocket follows 
a different pattern. In Fig. 3A, the F^ values are 
plotted against the tilt angle cj) of two SmBg single 
crystals. The uncertainty of F^ is determined by 
the half width at half-height of the F^ peak in the 
FPT plot. Similar to F^, the pattern has a 
fourfold symmetry, but the minima are located 
at ({) = 0° 90° and 180° that is, along the (100) 
crystalline axes. Fitting the data to the 2D form 
for this family, F^/cos((t)-0°), F^/cos((t)-90°), 
and F^/cos((t)-180°), with F^ = 30.5 T (solid 
lines in Fig. 3A) results in reasonable agreement. 
This suggests that the observed a pocket arises 
from a surface state on the (100) plane families. 
Similarly, Fig. 3C shows that the angular de- 
pendence of the y pocket F^ follows the func- 
tional form of Fj/cos((t)-0°), Fj/cos((t)-90°), 
and Fj/cos((t)-180°) with F^ -385 T, suggest- 
ing the two-dimensionality of the y pocket on 
the (100) surface plane. 

Further, given the small value and large un- 
certainly of F^, there is still a chance that an 
extremely elongated 3D ellipsoidal Fermi surface 
may fit the F^ versus (\> dependence. Experiments 
with cleaner SmBg crystals may resolve the issue. 

The angular dependence of the oscillation fre- 
quency suggests that the Fermi surface p is two- 
dimensional and likely arises from the crystalline 
(101) plane. In contrast, most of the theoretical 
modeling focuses on the surface states in the 




B 2.0 



0.02 0.04 



0.06 0.08 0.10 



0.12 0.14 



^ 1.6 



0 1.2 

■D 

1.0.8 
E 

CD 

H 0.4 



0.0 





(^ = 32° 




Y 






aJ 













200 400 600 800 

Frequency (T) 



1000 



Fig. 2. Quantum oscillation pattern observed by torque magnetometry. (A) The oscillatory mag- 
netic torque Xosc is plotted as a function of l/|ioH. Xosc is calculated by subtracting a polynomial back- 
ground from the full torque signal. (B) The Fast Fourier Transform (FFT) of the oscillatory torque Tqsc 
shows three major peaks, labeled a, (3, and y. The higher harmonics of the p peak are also observed, 
together with the p' peak from the neighboring surface plane. (C) The temperature dependence of the 
FFTamplitude of the three major peaks, normalized by theirO K limit. The magneticfield tilt angle^) ~ 32°. 
Fitting the oscillatory amplitude to the LK formula (solid lines) yields the effective mass m = 0.119me for 
the low-frequency a pocket, m = 0.129me for the high-frequency (3 pocket, and m = 0.192me for the 
higher-frequency y pocket. (D) At 0.3 K and <\) ~ 30°, the oscillation amplitude is tracked as a function of 
field H, generally known as the Dingle plot. Solid lines are fits to the Dingle factor. 
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(100) planes (31), in which the mapping of the 
band inversion X-point gives two Fermi pocket 
cross sections on the (100) surface plane. Our 
observed a and y pockets may be the two pre- 
dicted pockets. Recent angle resolved photo- 



emission spectroscopy (ARPES) measurements 
on SmBg revealed two Fermi pocket areas on the 
(100) surface Brillouin zone {32-36). Our mea- 
sured Fermi surface area and the mass of pocket 
y are comparable to those of a small Fermi sur- 
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Fig. 3. Angular dependence of the oscillating frequencies in SmBe. (A) Oscillation frequency 
of the observed Fermi surfaces a is shown as a function of 4). The error bar of P is defined as the half- 
width at half-height of the peak in the FFT plot. To improve the statistics, two SmBe samples were 
cooled down twice. All the oscillation frequencies are inferred from the FFTspectra at many tilt angles ^ 
from -30° to 180° (figs. S3 to S5) (37). The data are in agreement with a fit to the solid lines (described 
by Fo/cos((]) - 4)o). where = 0° 90° and 180°). The lines correspond to 2D Fermi surfaces on the 
surface planes (100), (001), and (-10 0). (B) Oscillation frequency F^of the observed Fermi surfaces 
p measured in two SmBe samples. is consistent with the expected 2D Fermi surface sitting at the 
four surface planes of the (101) surface family.The solid line is the fitting curve Fo/cos(4) - <\)o), where 
4)0 = -45°, 45°, 135°, and 225°. (C) The angular dependence of the oscillation frequency P of the 
pocket y. Solid lines are fits to Fo/cos(4) - (^o) with (^q = 0°, 90°, and 180°. (D) Sketch of the Brillouin 
zone of SmBe. Band inversion occurs at the X points. The projection of these X points to the (101) 
plane shows a possible location of the Fermi surfaces in the (101) surface state. 



Table 1. Parameters in the three Fermi pockets in SmBe- The oscillation frequency F and Fermi 
wavevector kf are obtained from the dHvA period. The effective mass m, Fermi velocity Vf, mean free 
path /, and the mobility \i are inferred from the Tand H dependences of the oscillation amplitude. The 
Landau level index plot yields the infinite field limit §. 





a 


P 


y 


F(T) 


29.3 ± 4 


286 + 5 


385 + 5 


Crystalline surface Origin 


(10 0) 


(10 1) 


(10 0) 


(nm"^) 


0.30 ± 0.02 


0.941 + 0.008 


1.080 + 0.011 


m/nrie 


0.119 ± 0.007 


0.129 + 0.004 


0.192 + 0.005 


Vf (10^ms~^) 


2.90 ± 0.39 


8.45 + 0.33 


6.50 + 0.21 


/(nm) 


22 + 8 


52 + 6 


25 + 4 


|i (xlO^ cmVV-s) 


1.1 + 0.5 


0.84 + 0.10 


0.36 + 0.05 


kfl 


6.6 + 2.9 


49 + 6 


27+4 


5 


-0.45 + 0.07 


-0.44 +0.06 


-0.32 + 0.07 



face centered at the r point and measured by 
ARPES. For a detailed comparison of our dHvA 
result and the ARPES results, see (37). 

Questions remain as to the origin of the ob- 
served Fermi surface in the (101) plane. In the 
(101) plane, there are four high symmetry points 
[(0, 0), (0, n/a), (n/V^a, 0), and (n/V^a, n/a)l 
For example (Fig. 3D), projecting the bulk band 
X points to the (101) plane leads to a pocket at (0, 
n/a) and a pair of pockets at (7r/\/2a,0) and 
(-71/ \/2a,0). Based on the time-reversal symme- 
try, only (0, n/a) has a Dirac point, v^hich is 
consistent with the experimental observation. 
However, the topological theory does not pro- 
hibit pairs of Dirac points at low symmetry points 
(37), which offers another possible origin for the 
observed pockets. 

The other important question is whether 
the 2D electronic state on the surface follows 
the Dirac dispersion. A general test is to track the 
Landau Level index plots to find out the infinite 
field limit, i.e., the geometric Berry phase factor. 
Using 45 T, the quantum limit is reached for 
pocket a, which, in the infinite magnetic field 
limit, points to -0.45 + 0.07, very close to -1/2, 
the geometric Berry phase contribution, simi- 
lar to other 2D Dirac electronic systems such as 
graphene (18, 19). As the oscillation frequencies 
of Fermi pocket p and y are quite close, filtering 
is needed to isolate the H dependence of the 
oscillation patterns for each Fermi pocket. This 
filtering may cause additional uncertainty of 
the oscillation phase, and the Zeeman effect and 
correlations may lead to some nonlinear effect 
of the Landau level index plots (37). 

Because Al flux is used in the sample growth, 
the observed effective masses, oscillation fre- 
quencies, and angular dependence bear some 
similarity to those of pure aluminum (Al) (38). 
However, the oscillation periods of all the ob- 
served Fermi surfaces are symmetric about crys- 
talline symmetry axes in the rotation plane. The 
observed fourfold symmetry and the behavior 
- l/cos(t), FP-l/cos((t)-45°) andi^ - l/cos(^ 
cannot be explained by a residual Al impurity 
[for more details, see (57)]. Therefore, the ob- 
served quantum oscillation pattern is an intrin- 
sic property of single-crystalline SmBg. 

We note also that our observed quantum os- 
cillation feature is quite robust against oxidiza- 
tion, as the samples were always in atmosphere 
for storage. The ordinary surface states knovm to 
occur for vacuum clean surfaces of metal hexa- 
borides such as LaBg disappear under even mod- 
est oxygen exposure (37, 39). 
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Pressure-dependent reactions are ubiqui- 
tous in applications of gas-phase chem- 
ical kinetics to practical problems, such 
as combustion (<100 bar), atmospheric 
chemistry (<1 bar), and chemical vapor 
deposition (<<1 bar). These reactions can take 
various forms involving chemical or thermal acti- 
vation, single or multiple unimolecular potential 
wells, and single or multiple sets of bimolecular 
products. Such reactions are enormously com- 
plicated problems to treat theoretically, even for 
a single-channel, single-well dissociation, which 
is the case we focus on here. The accurate first- 
principles prediction of pressure-dependent rate 
coefficients would dramatically improve the utility 
of theoretical kinetics as a tool for global chem- 
ical modeling. 

There are two constants of motion in an iso- 
lated dissociating molecule: the total internal 
(vibrational-rotational) energy, E, and the total 
angular momentum, J. A physically correct de- 
scription of the kinetics requires a quantitative 
understanding of the rate of collision-induced 
transitions in these two variables, kc{EJ',E\J'\ 
and the corresponding E- and J-resolved disso- 
ciation/isomerization rate constant, k{EJ). The 
master equation treats the time dependence of 
the state-resolved species populations arising from 
collisions and dissociation. This time dependence 
can be directly related to pressure-dependent 
thermal dissociation rate constants k{T,p) (1-4). 

Although many of the foundational formu- 
lations of reaction rate theory considered pressure 
dependence, subsequent theoretical developments 
have almost exclusively focused on pressure- 
independent reactions. Such pressure-independent 
rate coefficients are generally obtained from a 
thermal average of k(EJ) and are independent 
of the collisional transition process. Many classes 
of pressure-independent reactions can now be 

^Combustion Research Facility, Sandia National Laboratories, 
MS 9055, Livermore, OA 94551-0969, USA. ^Chemical 
Sciences and Engineering Division, Argonne National 
Laboratory, Argonne, IL 60439, USA. 
*Corresponding author. E-mail: sjk@anl.gov 



treated with very high accuracy. For example, 
quantitative first-principles kinetics calculations 
for pressure-independent bimolecular reactions 
of three- and four-atom (5-7) and larger systems 
(8) have been highlighted as milestones in the 
development of reaction rate theory. For larger 
systems, ab initio transition state theory provides 
high-accuracy predictions for the high-pressure 
addition rates as illustrated, for example, in our 
recent studies of 0(^P) + alkenes (9), of ^CHa + 
alkenes (10), and of alkyl radical + alkyl radical 
reactions (11), where a priori theory and experi- 
ment agreed to within 20%. In contrast, accurate 
first-principles methods for predicting pressure- 
dependent reaction rates have not been sim- 
ilarly developed. 

In the low-pressure limit, k(T,p) is indepen- 
dent of k(E,J) and is instead determined by the 
details of the collision-induced transitions in E 
and J (12). At intermediate pressures, both the 
collision-induced transition rates and the de- 
composition rates play an important role in the 
kinetics. Prior theoretical kinetics work has em- 
ployed simple empirical models for the colli- 
sional transition rates, with parameters adjusted 
to reproduce available experimental data (13, 14). 
This empiricism severely limits the utility of the- 
oretical kinetics as a mechanism development 
tool. Empirical theories are effectively limited 
to interpolations and extrapolations of the 
pressure and temperature ranges for reactions 
that have already been studied experimentally. 
Correlation of the energy transfer parameters 
with physical properties (e.g., molecular size, col- 
lision bath gas, reaction temperature) provides 
some predictive ability, but currently such cor- 
relations are not well determined and introduce 
additional uncertainties in the rate estimates of 
at least a factor of 2. For example, information 
regarding the temperature dependence of the 
energy transfer parameters is very limited, es- 
pecially for temperatures that deviate strongly 
from 300 K Furthermore, for complex multiple- 
well, multiple-product reactions, the effects of 
the uncertainties in such correlations are greatly 
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